Quantization Bounds on Grassmann Manifolds of 
Arbitrary Dimensions and MIMO Communications 
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Abstract — This paper considers the quantization problem on 
the Grassmann manifold with dimension n and p. The unique 
contribution is the derivation of a closed-form formula for the vol- 
ume of a metric ball in the Grassmann manifold when the radius 
is sufficiently small. This volume formula holds for Grassmann 
manifolds with arbitrary dimension n and p, while previous results 
are only valid for either p = 1 or a fixed p with asymptotically 
large n. Based on the volume formula, the Gilbert- Varshamov and 
Hamming bounds for sphere packings are obtained. Assuming a 
uniformly distributed source and a distortion metric based on 
the squared chordal distance, tight lower and upper bounds are 
established for the distortion rate tradeoff. Simulation results 
match the derived results. As an application of the derived 
quantization bounds, the information rate of a Multiple-Input 
Multiple-Output (MIMO) system with finite-rate channel-state 
feedback is accurately quantified for arbitrary finite number of 
antennas, while previous results are only valid for either Multiple- 
Input Single-Output (MISO) systems or those with asymptotically 
large number of transmit antennas but fixed number of receive 
antennas. 

I. Introduction 

The Grassmann manifold Gn,p (L) is the set of all p- 
dimensional planes (through the origin) of the n-dimensional 
Euclidean space L", where L is either R or C. It forms 
a compact Riemann manifold of real dimension f3p (n — p), 
where (3=1/2 when L = M/C respectively. The Grassmann 
manifold provides a useful analysis tool for multi-antenna com- 
munications (also known as Multiple-Input Multiple-Output 
(MIMO) communication systems. For non-coherent MIMO 
systems, sphere packings on the Gn,p (L) can be viewed as a 
generalization of spherical codes [l]-[3]. For MIMO systems 
with finite rate channel state feedback, the quantization of 
beamforming matrices is related to the quantization on the 
Grassmann manifold [4]-[6]. 

The basic quantization problems addressed in this paper 
are the sphere packing bounds and distortion rate tradeoff. A 
quantization is a mapping from the Gn,p (L) into a subset of 
the Q n:P (L), known as the code C. Define S = S (C) as the 
minimum distance between any two elements in C. The sphere 
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packing bounds relate the size of a code and a given minimum 
distance S. Assuming a randomly distributed source on the 
G n ,p (L) and a distortion metric, the distortion rate tradeoff is 
described by either the minimum expected distortion achievable 
for a given code size (distortion rate function) or the minimum 
code size required to achieve a particular expected distortion 
(rate distortion function). 

For the sake of applications [4]-[6], the projection Frobenius 
metric (i.e. chordal distance) is employed throughout the paper 
although the corresponding analysis is also applicable to the 
geodesic metric [3]. For any two planes P,Q £ Gn,p (IL), the 
principle angles and the chordal distance between P and Q are 
defined as follows. Let Ui e P and vi S Q be the unit vectors 



such that 



ujvi 



is maximal. Inductively, let U; e P and G 
Q be the unit vectors such that u\uj = and vjvj = for 



all 1 < j < i and 
defined as 8. 
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is maximal. The principle angles are 
for i =!,••• ,n [7], [8]. The 



chordal distance between P and Q is defined as 



dc{P,Q) = 
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The invariant measure on the Q n>p (L) is defined as fol- 
lows. Let 0(n)/U(n) be the group of n x n orthogo- 
nal/unitary matrices respectively. Let A S O (n) j\J (n) and 
B G O (n) /U (n) when L = R/C respectively. An invariant 
measure /i on the Q n ^ v (L) satisfies, for any measurable set 
■M. C G n ,p (L) and arbitrarily chosen A and B, 

fj. (AM) = \i (M) = n (MB) , 

The invariant measure defines the uniform distribution on the 

Gn,p (L) [7]. 

With a metric and a measure defined on the Q n ,p (L), there 
are several bounds well known for sphere packings. Let 5 be 
the minimum distance between any two elements of a code C 
and B (S) be the metric ball of radius 5 in the Gn,p (L). If K 
is any number such that K/j, (B (5)) < 1, then there exists a 



code C of size K + 1 and minimum distance <5. This principle 
is called as the Gilbert-Varshamov lower bound [3], i.e. 

1 
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(1) 



On the other hand, \C\fi(B (8/2)) < 1 for any code C. The 
Hamming upper bound captures this fact as [3] 

1 



\c\< 



(2) 



n(B(S/2)) 

These two bounds relate the code size and a given minimum 
distance 8. 

Distortion rate function gives another important property 
of quantization. Assume that Q is a random plane uniformly 
distributed on the Q np (L) and a distortion metric defined by 
the squared chordal distance d 2 . The average distortion of a 
given C is 



D (C) = E t 



Q 



min d' 

PGC 



(P,Q) 



(3) 



The distortion rate function gives the minimum average distor- 
tion for a given codebook size K, i.e. 



D* (K) = inf 

C:\C\=K 



D(C). 



(4) 



There are several papers addressing quantization problems in 
the Grassmann manifold. The exact volume formula for a B (8) 
in the Q np (C) where p = 1 is derived in [4]. An asymptotic 
volume formula for a B (8) in the Q n ^ p (L), where p > 1 
is fixed and n approaches infinity, is derived in [3]. Based 
on those volume formulas, the corresponding sphere packing 
bounds are developed in [3], [5]. Besides the sphere packing 
bounds, the rate distortion tradeoff is also treated in [9], where 
approximations to the distortion rate function are derived by the 
sphere packing bounds. However, the derived approximations 
are based on the volume formulas [3], [4] only valid for some 
special choices of n and p, i.e. either p = 1 or fixed p > 1 
with asymptotic large n. 

This paper derives quantization bounds for the Grassmann 
manifold with arbitrary n and p when the code size is large. An 
explicit volume formula for a metric ball in the Q n ^ p (L) is de- 
rived when the radius is sufficiently small. Based on the derived 
volume formula, the sphere packing bounds are obtained. The 
distortion rate tradeoff is also characterized by establishment 
of tight lower and upper bounds. Simulation results match the 
derived bounds. As an application of the derived quantization 
bounds, the information rate of a MIMO system with finite rate 
channel state feedback is accurately quantified for abitrary finite 
number of antennas for the first time, while previous results 
are only valid for either Multiple-Input Single-Output (MISO) 
systems or those with asymptotically large number of transmit 
antennas but fixed number of receive antennas. 

II. Metric Balls in the Q n , p (L) 

In this section, an explicit volume formula for a metric 
ball B (8) in the Q n<p (L) is derived. The volume formula is 
essential for the quantization bounds in Section [III] 



The volume calculation depends on the relationship between 
the measure and the metric defined on the Q n ^ p (L). For the 
invariant measure /i and the chordal distance d c , the volume of 
a metric ball B (8) can be calculated by 



(5) 
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where ,8 p are the principle angles and the differential 

form dfie is given in [7], [10]. 

The following theorem expresses the volume formula as an 
exponentiation of the radius S. 

Theorem 1: Let B (8) be a ball of radius 8 in Q n p (L). When 
8< 1, 



lx{B(8)) = 



_ / c n ^SP^Pl(l + o(8)) if L 



if L = 



(6) 



where (3 = 1/2 when L = M/C respectively and c„. P ( 3 is a 
constant determined by n, p and (3. When L = C, c„ iPi 2 can 
be explicitly calculated 
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When L = R, c n , P ,i is given by 
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and 



_ A 2 (p - i + 1) A (n - p - i + 1) 
W-ll 2A(n-i + l) 
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A(p) = 



The proof of Theorem Q] is not included due to the length 
limit. 

Theorem Q] provides an explicit volume approximation for 
real Grassmann manifolds and an exact volume formula for 
complex Grassmann manifolds when 8 < 1. Simulations show 
that this approximation remains good for relatively large 8 (Fig. 
®. 

Theorem Q] is consistent with the previous results in [4] and 
[3], which pertain to special choices of n and p and are stated 
as follows. 

Example 1: Consider the volume formula for a B (8) in the 
Qn.p (C) where p = 1. It has been shown in [4] that 

i i{B{8)) = 8 2 ^- 1 \ 
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(b) Complex Grassmann manifolds 
Fig. 1. The volume of a metric ball in the Grassmann manifold 

Theorem Q] is consistent with it where (3 = 2 and 0,1,1,2 = 1. 

Example 2: When p is fixed and n — > +00, the asymptotic 
volume formula for a B (5) is given by Barg [3] as 



r f3np+o(n) 



(9) 



On the other hand, Theorem Q] contains an asymptotic formula 
for L = C, 8 < 1, fixed p and asymptotically large n in the 
form 

/ § \ 2p("-p)+o(n) 



This follows from (0 and Stirling's approximation. Therefore, 
Theorem [T] is consistent with Barg's formula ([9). 

Importantly though, Theorem [T] is distinct from the previous 
results of [4] and [3] in that it holds for arbitrary n and p, 

1 < p < n. 

Fig. Q] compares the exact volume of a metric ball (0 and 
the volume evaluated by ©. For the volume approximation 
Cn,p,p$^ n ~ p \ the constant c„. P ) g is calculated either by (|7]i if 
L = C or by Monte Carlo numerical integral of (O if L = K. 
Simulations show that the volume approximation is close to the 
exact volume when the radius of the metric ball is not large. 



We also compare our approximation with Barg's approximation 
(8/^/pf np for n = 10 and p = 2 case. Simulations show 
that the exact volume and Barg's approximation may not be in 
the same order while the approximation in this paper is more 
accurate. 

III. Quantization Bounds 

Based on the volume formula given in TheoremQ] the sphere 
packing bounds are derived and the rate distortion tradeoff is 
characterized in this section. 

The Gilbert- Varshamov and Hamming bounds on the 
Q n>p (L) are given in the following corollary. 

Corollary 1: When 5 is sufficiently small, there exists a code 
in Gn,p (L) with size K and the minimum distance 8 such that 

n,p,/3° 

For any code with the minimum distance 8, 



K < 



*n,p,0 



-0p(n-p) 



Here and throughtout, the symbol < indicates that the inequal- 
ity holds up to (l + o (1)) error. 

Proof: The corollary follows by substituting the volume 
formula © into (HJ and (01. ■ 

The distortion rate function is characterized by establishing 
tight lower and upper bounds. 

Theorem 2: Let t = (3p(n — p) be the number of the real 
dimensions of the Grassmann manifold Q n:P {L). When K is 
sufficiently large, the distortion rate function is bounded by 



t 



(c n ,p,pK) ~<D*(K)< 



2T 



(c n ,p,f3K) 



(10) 



Due to the length limit, we only sketch the proof here. 
The lower bound is proved by an optimization argument. 
The key is to construct an ideal quantizer, which may not 
exist, to minimize the distortion. Suppose that there exists K 
metric balls of the same radius 80 covering the whole Q UtV (L) 
completely without any overlap. Then the quantizer which maps 
each of those balls into its center gives the minimum distortion 
among all quantizers. Of course such a covering may not exist, 
provding a lower bound of the distortion rate function. 

The upper bound is derived by characterizing the average 

distortion of the ensemble of random codes. Define a random 

code with size K as C ran d = {Pi, P2, ■ ' ' >Pk} where Pi's 

are independently drawn from the uniform distribution on the 

Q n . p (L). For any given Q <E Q n , P (L), define X, = d 2 c (Pi, Q) 

and W K = min (X x , ■■■ , X K ) = min dl (Pi, Q). Since 

Pi ec rand 

the codewords P^s 1 < i < K are independently drawn 
from the uniform distribution on the Q n ,p (L), Xi's 1 < i < 
K are independent and identically distributed (i.i.d.) random 
variables with the cumulative distribution function (CDF) given 
by Theorem Q] According to X t 's CDF, the CDF of W K can 
be calculated by extreme order statistics. We prove that for any 

2r (f) r -j 



given Q e G n , P (L), -E Wk [W k ] converges to — —^n, P ,0 
as K approaches infinity. Thus, K 2 • Eq \&w k \Wk\\ — ■ 
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Fig. 2. Bounds on the distortion rate function 



Ec ra „ d [D (Crand)] converges to the same constant, providing an 
upper bound of D* (K). 

It is worthy to point out that since the upper bound is 
corresponding to the ensemble of random codes, it is often used 
as an approximation to the distortion rate function in practice. 

Fig.|2]compares the simulated distortion rate function with its 
lower bound and upper bound in ( fTUb . To simulate the distortion 
rate function, we use the max-min criterion [5] to design codes 
and use the minimum distortion of the designed codes as the 
distortion rate function. Simulations show that the bounds in 
( [Tol l hold for large K. When K is relatively small, the formula 
( TTOb can serve as good approximations to the distortion rate 
function as well. In addition, we compare our bounds with 
the approximation (the "x" markers) derived in [9]. While the 
approximation in [9] works for the case that n = 10 and p = 2 
but doesn't work when n < 8 and p = 2, the bounds in ( TTOb 
hold for arbitrary n and p. 

IV. Application to MIMO Systems with Finite Rate 
Channel State Feedback 

As an application of the derived quantization bounds on 
the Grassmann manifold, this section discusses the information 
theoretical benefit of finite rate channel state feedback for 
MIMO systems using power on/off strategy. We will show that 
the benefit of the channel state feedback can be accurately char- 
acterized by the distortion of a quantization on the Grassmann 
manifold. 

The effect of finite rate feedback on MIMO systems using 
power on/off strategy has been widely studied. MIMO systems 
with only one on-beam are discussed in [4], [5], where the 
performance analysis is derived by geometric arguments in the 
Gn,i (C). For MIMO systems with multiple on-beams, many 
works, e.g. [9], [11], [12], employ Barg's formula (0 for 
performance analysis, which is only valid for MIMO systems 
with asymptotically large number of antennas but fixed number 
of receive antennas. Valid for arbitrary MIMO systems, the loss 
in information rate is quantified for high SNR region in [13], 
which is hard to be generalized to other SNR regions. For all 
SNR regimes, a formula to calculate the information rate is 
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Fig. 3. System model 



proposed in [14] by letting the numbers of transmit and receive 
antennas and feedback rate approach infinity simultaneously. 
But this formula overestimates the performance in general. 

The system model of a wireless communication system with 
Lt transmit antennas, Lr receive antennas and finite rate 
channel state feedback is given in Fig. [3] The information bit 
stream is encoded into the Gaussian signal vector X € C sxl 
and then multiplied by the beamforming matrix P 6 C LtXS 
to generate the transmitted signal T = PX, where s is the 
dimension of the signal X satisfying 1 < s < Lt and the 
beamforming matrix P satisfies Ptp = I s . In power on/off 
strategy, E [XX*] = P on I s where P on is a positive constant to 
denote the on-power. Assume that the channel H is Rayleigh 
flat fading, i.e., the entries of H are independent and identi- 
cally distributed (i.i.d.) circularly symmetric complex Gaussian 
variables with zero mean and unit variance (CAf (0, 1)) and H 
is i.i.d. for each channel use. Let Y € C LflXl be the received 
signal and W e t£ L R xl be the Gaussian noise, then 

Y = HPX + W, 



where E [WW+] = I Ln . We also assume that there is a 
beamforming codebook B = {P; S C LtXS : PjP; = I s | 
declared to both the transmitter and the receiver before the 
transmission. At the beginning of each channel use, the channel 
state H is perfectly estimated at the receiver. A message, which 
is a function of the channel state, is sent back to the transmitter 
through a feedback channel. The feedback is error-free and rate 
limited. According to the channel state feedback, the transmitter 
chooses an appropriate beamforming matrix P^ G B. Let the 
feedback rate be i?fbbits/channel use. Then the size of the 
beamforming codebook \B\ < 2 Rtb . The feedback function is 
a mapping from the set of channel state into the beamforming 
matrix index set, cp : {H} — * {i : 1 < i < \B\}. This section 
will quantify the corresponding information rate 



X = max max E 

B:|e|<2 ii fb V 



log 



Pc 



' HP WH) P ^(H) H 



where P on = p/s and p is the average received SNR. 

Before discussing the finite rate feedback case, we consider 
the case that the transmitter has full knowledge of the channel 
state H. In this setting, the optimal beamforming matrix is 
given by P opt = V s where V s € C il-Xs is the matrix 
composed by the right singular vectors of H corresponding to 
the largest s singular values [6]. The corresponding information 
rate is 



2o P t = Eh 



^ln(l + P on Ai) 



(11) 
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Fig. 4. Performance of finite size beamforming codebook. 

where A,; is the i th largest eigenvalue of HHt In [6], we derive 
an asymptotic formula to approximate a quantity of the form 
Eh Ei=i m (1 + c ^i)] where c > is a constant. Apply the 
asymptotic formula in [6]. T op t can be well approximated. 

The effect of finite rate feedback can be characterized by the 
quantization bounds in the Grassmann manifold. For finite rate 
feedback, we define a suboptimal feedback function 

i = <P (H) = arg min d\ (V (Pi) , V (V,)) , (12) 

l<i<|B| 

where V (P^) and V (V s ) are the planes in the Gl t .s (C) 
generated by Pi and V s respectively. In [6], we show that this 
feedback function is asymptotically optimal as i?fb — > +oo and 
near optimal when Rib < +oo. With this feedback function and 
assuming that the feedback rate Rfo is large, it has been shown 
in [6] that 



E 



H 



^ln(l +?7sup-PonAi 



i=l 



(13) 



where 



= 1 - - inf E Vs 

SB:|B|<2 R fb 

= 1 - -D* (2 Rlb ) . 



min di {V (Pi),V(y B )) 

1< % < LB 



(14) 



Thus, the difference between perfect beamforming case (fTTT i 
and finite rate feedback case (Tl3T > is quantified by ?7 SU p, which 
depends on the distortion rate function on the Gl t ,s (C). Substi- 
tute quantization bounds ( fT0T > into ( TT4T > and apply the asymptotic 
formula in [6] for Eh Ei=i m (1 + c ^i)]- Approximations to 
the information rate X are derived as functions of the feedback 
rate R^. 

Simulations verify the above approximations. Let m = 
min (Lt, Lr)- Fig. |4] compares the simulated information rate 
(circles) and approximations as functions of R^jm 2 . The infor- 
mation rate approximated by the lower bound (solid lines) and 
the upper bound (dotted lines) in (TTOb are presented. As a com- 
parison, we also include another performance approximation 
(dash-dot lines) proposed in [14], which is based on asymptotic 



analysis and Gaussian approximation. The simulation results 
show that the performances approximated by the bounds ( TTOb 
match the actual performance almost perfectly and are much 
more accurate than the one in [14]. 

V. Conclusion 

This paper considers the quantization problem on the Grass- 
mann manifold. Based on the explicit volume formula for 
a metric ball in the Q n ^ p (L), the corresponding Gilbert - 
Varshamov and Hamming bounds are obtained. Assuming 
the uniform source distribution and the distortion defined by 
the squared chordal distance, the distortion rate function is 
characterized by establishing tight lower and upper bounds. As 
an application of these results, the information rate of a MIMO 
system with finite rate channel state feedback is accurately 
quantified for abitrary finite number of antennas for the first 
time. 
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